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ABSTRACT 

We examine the thermal free-free mihimetre fluxes expected from non-dusty and non- 
pulsating K through mid-M giant stars based on our limited understanding of their 
inhomogeneous chromospheres. We present a semi-analytic model that provides es- 
timates of the radio fluxes for the mm wavelengths (e.g., CARMA, ALMA, JVLA 
Q-band) based on knowledge of the effective temperatures, angular diameters and 
chromospheric Mg II h & k emission fluxes. At 250 GHz, the chromospheric optical 
depths are expected to be significantly less than unity, which means that fluxes across 
the mm and sub-mm range will have a contribution from the chromospheric mate- 
rial that gives rise to the ultraviolet emission spectrum, as well as the cool molecular 
material known to exist above the photosphere. We predict a lower bound to the in- 
ferred brightness temperature of red giants based on heating at the basal-flux limit 
if the upper chromospheres have filling factor ~ 1. Multi- frequency mm observations 
should provide important new information on the structuring of the inhomogeneous 
chromospheres, including the boundary layer, and allow tests of competing theoretical 
models for atmospheric heating. We comment on the suitability of these stars as mm 
flux calibrators. 

Key words: Stars: chromospheres - Starsdate-type - Stars: activity - Radio: con- 
tinuum stars. 



1 INTRODUCTION 

Convective motions beneath the photospheres of late-type 
(cool) stars agitate the plasma, leading to the excitation 
of acoustic and magnetic disturbances. These waves induce 
non-radiative heating of the upper photosphere and chro- 
mosphere. Essentially all cool stars possess material heated 
above the prediction of radiative equilibrium, which emits 
a chromospheric spectrum that is most conspicuous in the 
ultraviolet (uv). The specific roles of acoustic and magnetic 
waves in heating the quiet solar chromosphere and the role of 
magnetic fields in structuring its atmosphere are the subject 
of considerable de bate ([Kalkofen, Ulmschneider, fc Avrett| 



1999 Ayres 2002 Vecchio, Cauzzi, & Reardon 20091. For 



the Sun and red giants the chromospheric plasma has a 
detectable signature in the millimetre (mm) thermal con- 



tinuum (Altenhoff, Thum, & Wendker 1994 Loukitcheva, 
Solanki, Carlsson, fc Stein|2004[ ) 



Stellar chromospheric surface fiuxes measured from 
dominant emission lines, such as Ca II H & K and Mg II 
h & k, show a wide scatter at a given effective tempera- 
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ture (Tefl), just as found for regions of different activity on 
the Sun. However, there is a well defined empirical lower 
bound known as the basal-flux, which is a strong function of 
effective temperature, -Fbasai ~ T^s, but is only a weak func- 
tion of surface gravity {g, ) . The chromospheric cooling from 
the quietest solar regions, in the centres of supergranules, is 
close to the basal-flux for a Taft = 5770 K star (e.g., Schri- 



jver & Zwaan 2000). Stars that have chromospheric emission 



above the basal-flux are thought to have additional contri- 
butions from areas of enhanced surface magnetic activity, 
analogous to solar active regions (sunspots and associated 
plage). However, no consensus has been reached about the 
origin of the basal-flux. One line of thought is that it arises 
from the deposition of purely acoustic shock energy gener- 



ated in the sub-photospheric convection zone ( Ulmschneider 



1991 1, while the absence of appropriate spectroscopic shock 



signatures in inactive basal-flux red giants ( Judge fc Carpen- 



|ter,1998 ) suggests a magnetic origin. In fact, on the Sun both 
mechanisms are seen to operate: small-scale transient shock 
heating in the so-called "K grains" and more steady mag- 
netic heating in the large-scale supergranulation network. 

The physical structures of chromospheres from these 
two mechanisms are likely to be quite different. The acous- 
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tic wave picture would lead to a 3-D network of shocks that 
is highly time variable and at a given position the gas tem- 
perature would fluctuate from very cool to very hot with 



20081, whereas studies of inactive basal-flux red giants may 



the mean temperature being cool ( Wedemeyer-Bohm et al 



20071. This mechanism leads to localized intermittent chro- 



mospheric emission but, which when averaged over the stel- 
lar disk, would not be apparent in the mm-radio. The pres- 
ence of magnetic fields can lead to longer lived structures 
and enhanced atmospheric heating. Time-independent 1-D 
semi-empirical models for the Sun that represent regions of 
different activity are well established fVernaz za et al.|1976 



1981 Fontenla et al.|1990 |. Similar but spatially unresolved 
models have been constructed for a few cool stars includ- 
ing inactive red giants. These semi-empirical models are de- 
signed to reproduce the temporally and spatially averaged 
uv and optical emission line profiles and fluxes. A charac- 
teristic of these models is that they require a gradual tem- 
perature rise from the top of the photosphere to the up- 
per chromosphere, and then a very steep rise through the 



chromosphere-coronae transition region (Kelch et al. 1978 



Harper|[T992l [Luttermoser et al.|1994| |McMurry||1999 

One might expect that acoustic models, whose mean 
temperatures do not increase with height, would fail to pre- 
dict the same uv fluxes as semi-empirical models with an 
outward temperature increase, but this is not necessarily 
the case. In the uv, hv/kT 3> 1 and the source functions 
and emissivity oc e"''"'"^^ are very sensitive to tempera- 
ture. In this situation, the shock peaks can dominate the 
temporally and spatially averaged emission of the dynamic 
atmosphere. Thus it is possible for acoustic shock and semi- 
empirical models to produce the same chromospheric fluxes 
even though the temperature structures are profoundly dif- 
ferent ( jCarlsson fc Stein|1995 |. 

However, in the mm-radio, the thermal continuum 
source (Planck) function depends linearly on Tc, so both 
solar and stellar mm observations can probe differences be- 
tween intermittent and long-term structured atmospheres 
and potentially differentiate between acoustic and magnet- 
ically heating models. A comparison between uv and radio 
diagnostics could test the symbolic inequality 



-hv/kT 



-hv/k(T) 



(1) 



where the averages are over space and time. The LHS ex- 
emplifies uv emission diagnostics while the RHS reflects the 
same source function but with the mean temperature, (T), 
as inferred for example from mm radio emission. In prac- 
tice uv and mm-radio fluxes calculated from different atmo- 
spheric models can be compared to observations. The hy- 
drogen ionization determines the electron densities which, 
in part, determine both the mm-radio opacity and the uv 
collisional excitation rates. The remaining difference in the 
temperature dependence for formation leads to the test of 
equation ([TJ. 

Here we analyze published thermal mm (250 GHz) 
fluxes | |Altenhoff et al.|1994 1 from single non-dusty and non- 
pulsating K and M giants (luminosity class HI) to investi- 
gate their chromospheric signatures and explore the role of 
future multi-frequency observations in revealing the nature 
of chromospheric structuring. Future solar Atacama Large 
Millimeter/submillimeter Array (ALMA)imaging will be fo- 
cussed on obtaining sufficient spatial resolution to separate 



allow us to avoid the contamination from magnetic active 
regions by observing giants that are covered with inactive 
chromospheric regions analogous to quiet solar supergranule 
centres. 

A first estimate for the expected mm flux from a red gi- 
ant is provided by the hard disk model; which is appropriate 
when the ionized density scale height, //p°" is much smaller 
than the stellar radius, The observed flux density {Fi,) 
is given by 



(2) 



where F* is the average flux density at the surface, having 
angular diameter </> with respect to the distance observer. 
Working in terms of brightness temperatures Tbi we have 



Fv = 



2c2 



Tbi (m) ^J.d^^ T^kv (f) ^^^^^ 



(3) 



where is the usual cosine of the angle between a ray and 
the normal to the emitting surface, and (Tbt) is the weighted 
mean brightness temperature. In more convenient units, the 
flux in mJ^l^is 



Fv (mJy) = 13. i 



3900 K 



V 

250 GHz 



V 10 mas/ ^ ' 



Many nearby red giants now have measured angular di- 
ameters, (/>«, and estimated effective temperatures, Tofi 
( [Mozurkewich et al. 2003 1. Inserting 0* and substituting Tefi 
for (TBr) into equation ([4|) shows that many of these giants 
can in principle be observed at multiple frequencies with 
high signal-to-noise ratios when ALMA is completed given 
the expected sensitivity limits (Butler & Wooten 19991. 



Here we develop a more physically motivated semi-analytical 
model for mm fluxes that incorporates constraints on chro- 
mospheric structure from independent IR and uv observa- 
tions, which also allows us to include the effects of different 
levels of chromospheric heating. 

In §2 we compile pertinent empirical constraints derived 
from semi-empirical models and other studies of red giants; 
and in §3 estimate the expected chromospheric mm optical 
depths. In §4 we develop a semi-analytic model for {Tsr) 
based on these empirical constraints, which we then cali- 
brate against direct measurements of the well studied red 
giant a Tau (K5 III) (Aldebaran). A comparison with ob- 
servations and speciflc predictions for ALMA frequencies are 
presented in §5, and a discussion of the diagnostic potential 
of mm observations in 56. Conclusions are drawn in 57. 



2 EMPIRICAL CONSTRAINTS ON STELLAR 
CHROMOSPHERES 

In this section we gather together nuggets of information 
from previous studies of red giant chromospheres, which we 
will incorporate into our analytic model for the mm radio 
fluxes. The radio optical depth, r^, is proportional to differ- 
ent powers of the atmospheric properties - electron density, 
temperature, and path length - and we seek to estimate these 



the magnetic and non-magnetic regions (Loukitcheva et al 



1 mjy= lO-^^W m-'^Hz- 



10" 



^Hz- 
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for stars covering a range of T^s and different levels of stellar 
activity. 

Historically, the usual lack of temporal and spatial res- 
olution in observations of the chromospheres of distant stars 
partly motivated the development of 1-D time-independent 
semi-empirical models. The earliest models used diagnostics 
observable from the ground - Ca II H & K, the Ca II IR 
triplet, and Hq - and early satellite observations of Mg II 
h & k from Copernicus ( Ayres et al.||1974 Ayres & Linsky 



1975[ ). These types of models are analogous to the venera- 
ble VAL series of solar models, although for the latter spa- 
tial resolution did permit construction of components rep- 
resentative of regions of different activity levels ( |Vernazza| 



et al. 1976 1981). More recent variants of these solar com- 



ponent models are the FAL series (Fontenla et al. 19901, 
which are employed in solar irradiance studies ( |Fontenla| 
et al.|2009| l, among other uses. Two notable features of these 
semi-empirical models are: (1) the temperature inversion 
in the low chromopshere above a temperature minimum, 
T^^^ ~ O.TSTeff; and (2) that the electron density, ria, is, 
to within a factor of ~ 2, approximately constant through- 



out the chromosphere (Ayres 1979 Harper 19921. At the 
temperature minimum, the electron density (~ 10~''nH) is 
dominated by (photo) ionized metals; but as the tempera- 
ture increases through the chromosphere, hydrogen gradu- 
ally becomes more ionized until at the top of the chromo- 
sphere (~ 10* K), hydrogen is predominantly ionized and 
then no ~ nn. 

Thus, through the chromosphere, the increasing ion- 
ization of hydrogen counteracts the rapid outward decline 
of the density due to hydrostatic equilibrium, maintaining 
a more-or-less constant rie with increasing altitude (Ayres 



1979). However, once the hydrogen has become mostly ion- 
ized, there no longer is a ready source of bound electrons to 
be liberated, and so in the higher layers the electron density 
will fall off rapidly, following the hydrostatic decline in tih 
itself ( |Ayres|1979t . 

While there is a debate on whether the Sun has a per- 
manent chromospheric temperature rise or not - it should 
be noted that presently the semi-empirical models neverthe- 
less are the best thermodynamic inventory of stellar chromo- 
spheric plasma available. They approximately reproduce the 
amounts of plasma at different electron temperatures. This 
is discussed further in §6. 

Pre-IUE, Mg II h & k observations revealed a strong de- 
pendence of chro mospheric heating with effective tempera- 



ture, F mkii ~ T°s (Linsky & Ayres 1978 Ayres|l979l. Ayres 



(1979) developed a scaling law for the mass column density 



at the onset of the chromospheric temperature rise, which 
can be rephrased in terms of the chromospheric electron den- 
sity, rio, where atmospheric heating is sufficient to overcome 
H~ cooling, namely 



, /ll/2Trl/2„l/2^5/2 
■^Fc 9* -'eft • 



(5) 



^Fo represents the abundance of low first ionization poten- 
tial elements relative to solar, J" is the possible enhancement 
of a particular star relative to the general Tea scaling related 
to differences in intrinsic activity. In Ayres (1979), J- = 1 
was taken to represent levels appropriate to inactive stars 
or the quiet sun, and = 10 for solar plage regions or (fast 
rotating) stars showing enhanced activity. Some empirical 
evidence for this trend on ric in red giants is given by |Byrne| 



et al.| ( |1988| based on lUE C II] 23 25 A intersystem flux 
ratios, and by Buchholz et al.| ( |1998[ ) who compared time- 
dependent acoustically heated atmospheric models with the 
mass column density of an effective T^^. Chromospheric 
mm radio fluxes are sensitive to the opacity, which in turn 
depends on the square of the electron density. Here we as- 
sume the scaling in equation ([sj provides an appropriate 
dependence of ric on star specific parameters. 



2.1 Basal Fluxes 

We begin by considering the ratio of a star's measured Mg II 
h & k flux to the basal-flux, i.e., a basal-flux star has J- = 1. 
We recall the usual implicit assumption that Mg II h & k line 
fluxes represent a fixed fraction of the total chromospheric 



radiative losses (Linsky & Ayres 1978). The other major 



coolants are H Lya, Ca II, and Fe II (e.g., Judge & Stencel 
|1991| ). 

For this work we take lUE Mg II h & k fluxes from 



Martinez et al. (2011), which have typical uncertainties of 
±10%, and we also adopt their V — K -based expression for 
the basal- flux (erg cm'^ s""*^), 

log J'basai = 6.78 log Tcft - 19.74 (6) 



which has a very similar Tcft-dependence to [Linsky fc Ayres] 
( 1978| ) . We have chosen specifically their V — K form of the 
basal-flux, because the implied angular diameters are more 
reliable than the B — V form at lower TeS ■ [Judge fc Stencel] 
( |1991[ ) presented Mg II h fc k surfaces fluxes for a sample 
of low and intermediate mass giants, and their lower-bound 
drops below equation ||6| for TcS < 3300 K. 



2.2 Inhomogeneous Chromospheres 

Observations of the 4.6/im CO fundamental bands have 
demonstrated that on the Sun and late-type stars there 
exist regions of cool molecular gas with temperatures well 
below the T^y„ derived from ultraviolet and optical diag- 
nostics ( Ayres fc Testerman[[198T Wiedemann et al.|1994 ). 
On the sun the area fllling factor of the CO material has 
been proposed to be up to 80% in the low chromosphere 
( [Ayres ^2002) , but that factor is controversial. For the red 
giants the fllling factor is close to unity with the tempera- 
ture scaling as Too ~ 0.60roff (which is close to radiative 
equilibrium values when molecular cooling is included). At 
present there is a vigorous debate as to the exact nature 
of the CO clouds given the highly dynamic, time-dependent 



nature of the solar atmosphere ( 


Kalkofen, Ulmschneider, fc 


Avrett|1999| [Ayres[[2002[ [Leenaarts, Carlsson, Hansteen, fc 


Gudiksen[[2011 


). The presence of the CO material must be 



accounted for in chromospheric models, particularly when 
these regions are probed by mm continuum emission. For 
example, recent spectral interferometry of the much later 
spectral-type M giant BK Vir (M7 HI) has revealed the 
presence of spatially extended CO with an inner extended 
component at R — 1.27?,, with Tco — O.OSTofi ( [Ohnaka] 
et al.|2012| ). 



At lower effective temperatures, other molecules, such 
as water vapour ( Tsuji[ [2008), can play a role in further re- 



ducing the effective boundary layer temperature. [Quirren^ 



bach et al. ( 1993 ) found evidence from the Mark HI inter- 
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ferometer that the TiO continuum becomes significantly ge- 
ometrically extended by M5 III (R Lyr), for which the ratio 
of TiO to continuum diameters is 1.12±0.02. These extended 
molecular components possibly are related to the 'molsphere 
phenomenon' that is more common in the Mira and lower 
gravity M stars ( |Tsuji et al.|1997[ |Perrin et al.||2004 l. 



2.3 Chromospheric Extent 

The apparent angular size of a stellar radio source clearly is 
an important determinant of the radio flux. Here we consider 
the factors that control the angular diameter in the mm 
band, and in particular its ratio, G, to the photospheric 
angular diameter, i.e., (/>mm = G(j)t. 

The ratio of pressure scale- height, i/pros, to stellar ra- 
dius, -R, , is 



(7) 



so that red giants with their much larger radii, but near solar 
masses, have relatively thicker chromospheres than the Sun. 

Eclipse observations of the C, Aurigae systems, which 
have intermediate mass K lb primaries, show large plasma 
turbulence and direct evidence for atmospheric extensions 



greater than the thermal pressure scale height, (e.g., Eaton 



1993). This indicates that iifpics is controlled by both ther- 



mal and non-thermal motions. Here we describe the increase 
in -ffpros from non-thermal flows by a quasi-hydrostatic equi- 
librium, 



Hp 



2(/i)mH%* 



(8) 



where u is a characteristic velocity. For C, Aur systems, v 
is consistent with the turbulence inferred from spectra, and 
possibly correlated on spatial scales comparable to the den- 
sity/pressure scale-height, (p) is the chromospheric average 
of the mean mass per particle in units of the hydrogen mass 
ma, and is the surface gravity. We adopt (/i) — 1.2 corre- 
sponding to a partially ionized hydrogen [uc — Clnn) and 
surface helium abundance of whe = O.Itih. Semi-empirical 
models typically include this non-thermal term in the equa- 
tions of quasi-hydrostatic equilibrium. For Aur systems, 
the implied turbulence is slightly greater than the hydrogen 
sound speed, and in our model we set v = \/2kT. While esti- 
mates of stellar radii are available from interferometry (e.g., 
Baines et al.|20l"0 1 and nearby red giants also have Hippar- 
cos parallaxes, the stellar masses remain uncertain for the 
single stars. Typical estimates are 1 — 2Mq , leading to some 
uncertainty in the surface gravity in equation ([8|. 

The semi-empirical one-component model of a Tau 
| McMurry|1999[ ) predicts a 250 GHz (1.2 mm) radial optical 
depth of unity at 7? ~ 1.037?*, but the apparent angular di- 
ameter (half central intensity) corresponds to 0mm — 1.10*. 
This is a result of the larger column density of the tangen- 
tial sight-lines that deflne the 1.2 mm limb. The apparent 
angular diameter of a Boo predicted from a semi-empirical 
model of (Drake 19851 is i^mm — 1.06(/!>*, and for g Her 
(M6 HI) 4>mm ^ 1.08<^* ([Luttermoser et al.||T994|). These 



models were constructed with different diagnostics and as- 
sumptions about the contribution of turbulent pressure. So 
although there is a trend of increasing fractional chromo- 
spheric extension with later spectral type, the lower Tcfi and 



thus lower chromospheric heating, leads to lower and re- 
duced mm-optical thickness. In light of these competing fac- 
tors we adopt a typical extension of LOScp* (i.e., G = 1.08), 
and accept an additional uncertainty of ~ 2% in G which is 
consistent with the level of approximations adopted here. In 
summary, here (Tbi) for the mm radio is defined as 



2c^R 



7rfci/2 (G<^.)' 



(9) 



Very recently, single baseline, spectrally-resolved, K gi- 
ant chromospheric Ca II IR-triplet visibilities from CHARA 



have been presented by Berio et al. (2011 1. The authors at- 



tempted to estimate the chromospheric extent by using a 



correction procedure based on the Eriksson et al. ( 1983 1 



semi-empirical chromosphere of the coronal star /3 Get 



(KO HI). Berio et al. (2011 1 derive chromospheric extents of 
1.16 — 1.477?*, although these values are based on a model 
whose actual spatial extent is only 1.027?*. The extended 
emission observed in the cores of the IR-triplet lines depends 
on the presence of Ga II ions which might not be related to 
the regions of higher electron density to which the radio 
is sensitive. The calculation of the IR-triplet visibilities re- 
quires a careful treatment of the Ga II to Ca III photoioniza- 
tion by H Lya ( [Linsky fc Avrett||1970} |Rowe||1992[ ), spher- 
ical geometry, and possible impact of cross-redistribution 
of photons scattering within the Ga II H & K resonance 
lines (Uitenbrook 1989). The new GHARA measurements 



are very important diagnostics and one of us (NO'R) is un- 
dertaking these non-LTE radiative transfer calculations at 
the present time. 



3 CHROMOSPHERIC MILLIMETRE OPTICAL 
THICKNESS 

As mentioned earlier the chromospheric electron density has 
contributions from the low first ionization potential elements 
and from hydrogen. At low temperatures (T < 6000 K), ion- 
ization of hydrogen is a two stage process: excitation to the 
n = 2 level by electron collisions or scattered H Lya pho- 
tons followed by photoionization by the optically thin pho- 
tospheric Balmer continuum. At higher temperatures, direct 
collisional ionization becomes important. When hydrogen is 
partially ionized, free- free (thermal bremsstrahlung) opacity 
dominates at mm wavelengths. To examine the influence of 
the stellar outer atmosphere on the mm thermal emission, 
we first consider the optical thickness of the warm chromo- 
spheric component. For typical conditions, the opacity cor- 
rected for stimulated emission is given by |Rybicki fc Light^ 
manl (|2004|) 



0.01772 n^UionZC, 



-Qff cm 



(10) 



where the variables are in cgs units, Zion is the charge of 
the ions, ric and nion are the number densities of the elec- 
trons and ions, respectively, and ga is the free- free Gaunt 
factor. At cm- wavelengths the familiar power-law approxi- 



mation from Altenhoff et al. ( 1960 1 is 



cm -I -I —0.1 

ga = 11.96Te u 



(11) 
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Figure 1. Schematic diagram of an idealized static two-component inhomogeneous red giant chromosphere. The warm chromosphere 
temperature rises outwards from Tint to Tchrom E'nd had a total optical thickness Ar. Beneath the warm component, lies an opaque cool 
molecular region with Tqq < T^^^. 
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Figure 2. 250 GHz brightness temperatures for the sample given in Table 1. We also show three additional stars (fi Cep, t Aur, and 
7 Aql) which are of higher luminosity and might have emission in excess to that predicted by our model. The basal (Tar) predicted from 
our model is also shown as a solid curve. 



but at ALMA frequencies, which range from 84-950 GHz 
(Band 3-10), the power-law dependence is slightly different, 
namely 



sub — mm 

as 



= 24.10T+''-^'^i^" 



(12) 



as derived using the accurate Gaunt factors from |Hummer| 
(19881. Under typical chromospheric conditions, the major- 
ity of abundant species are either neutral or single ionized 
so that Zioii = 1 and Tie = Wion. This leads to an expression 
for the mm opacity. 
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0.427nS 



(13) 



The optical thickness of a slab of depth, AL, is then simply 
Atv = KvAL. 

The assumption that the electron density is approxi- 
mately constant, as the ionization nc/nu goes from 10~* — > 
1, implies that the hydrogen density has declined by four 
orders of magnitude and therefore the physical thickness of 
the layer AL ~ 9 density scale-heights. To the level of ap- 
proximation here we set Hp ~ Hprcs 

For the radial optical depth of the chromosphere we 

have 



At, 



0.427 < 



9- 



2kn 
{fj,)mHg, 



6.4 < k 



(14) 



This expression reveals only a very weak dependence on 
temperature variations within the chromosphere which is a 
result of the cancellation resulting from the combination of 
opacity and scale-height. In the following, we set the surviv- 
ing temperature factor to T^f^'^- Inserting into equation ( 141 



typical values for the M2.5 III star /3 Peg pecin et al.|2603 
M. = 2.OM0, R = llOi?*, Tcfi = 3600 K) and a ty pical uv 



based value of Jie ~ 6 x 10* cm ^ ( Byrne et al. 1988 1 , we find 



for 1.2 mm (250 GHz) and 3 mm (100 GHz) chromospheric 
optical thicknesses of At25oghz ~ 0.1 and Atiooghz — 0.8, 
respectively. Thus depending on the level of stellar activ- 
ity, the mm wavelength regime is one where chromospheres 
can go from optically thick at mm and cm wavelengths to 
optically thin in the sub-mm. 

If the optical thickness of the chromosphere is less than 
unity the next emitting layer we must consider is the cool 
molecular gas immediately above photosphere with its near 
unity filling factor. In these cool, high density layers, where 
the hydrogen density increases rapidly with depth, H~ free- 



free can become an important opacity source ( Harper et al 



20011. Here we assume that the bifurcated molecular re- 
gions are opaque at these wavelengths. Thus the optical 
depth unity is reached when Te ~ Tco, which is reason- 
able if the temperatures in these outer radiative equilibrium 
layers change only slowly with optical depth. 



4 DESCRIPTION OF THE SEMI-ANALYTIC 
MODEL 

We now assemble the different aspects discussed above and 
construct a semi-analytic model to provide estimates of 
the mm, e.g., ALMA band, radio fluxes. We then calibrate 
against well studied a Tau (K5 III) to normalize all the 
approximation constants of order unity accumulated along 
the way. A schematic representation of the model inhomo- 
geneous chromosphere is shown in Figure 1. Above the uv 
temperature minimum TmIii li^s the opaque cool molecu- 
lar material with temperature Too • Immediately above this, 
perhaps in an overlying magnetic canopy, is the base of the 
warm chromosphere with temperature Tint; Te then rises 
linearly in mass column density to the upper chromosphere 
with temperature Tchrom and has an optical thickness At. 
In the model is approximately constant (to within a fac- 
tor of ±2) and Tc, in the warm component, is approximately 
linear with height. The source function 5*^ oc Tc, and since 
the mm optical depth is oc ri^AL, then 5*^ is a linear 



function of . In this case the integration of the formal so- 
lution of the radiative transfer equation for the brightness 
temperature can be done analytically. 

Tbt, a function of frequency and viewing angle, repre- 
sents a weighted contribution of the three temperature pa- 
rameters and is given by 



(15) 



where Tco, Tint, Tchrom are the key temperatures described 
above and weighting coefficients are obtained from the for- 



mal solution of the transfer equation ( Olson & Kunasz 1987 1 , 
which in our notation are 



/? = 1 - (1 - q) /At 

7 = 1 — Q — /3. 



(16) 
(17) 
(18) 



From Wiedemann et al. ( 1994 \ we adopt Tco — O-dT^a 



and Tint = 4500 K. We take Tchrom^^ 10000 K where hy- 
drogen has become 50% ionized ( jMcMurry 1999) and any 
further increase in temperature now leads to a rapid decline 
in rie and hence also k,, creating a transparent boundary 
above Tchro m • 

In this analytic model, the radio fluxes are dependent 
on Teft and the optical depth of the warm chromosphere 
layer above the CO zone. If the emitting layers are not too 
thick compared to the stellar radius, the stellar flux can be 
described as the integral of Tsr (m) for a plane-parallel slab, 
where /i — cos 6 and 6 is the angle between the atmosphere 
normal and the ray. This allows us to write a weighted mean 
brightness temperature, e.g., in equation Q, as 



(TBr) 



Tbi (m) M dfi 



2 Tbi im) yawi = TBr (m = 1/2) 



(19) 



In keeping with the simple nature of our analytic model we 
replaced the integral by a single point Gaussian quadrature 
(n = 1, /ii = 1/2, and lui — 1). (Tbi) is found by substitut- 
ing At from equation 1 14 1, adjusted for the viewing angle 
p = 1/2, into equations (I61, (17l, and (I8I to evaluate the 
weighting coefficients and then solve for equation (15 1. 



4.1 Calibrating with a Tau (K5 III) 

To mitigate some of the inherent uncertainties in the model 
we evaluate the chromospheric optical thickness (equa- 
tion 14 1 for a Tau using empirical constraints. We initially 



use the scheme outlined above to calculate the 250 GHz flux 
and compare it to the measured 1986 IRAM [Altenhoff et al.| 
(1994) ) flux of 51 ± 6mJy. At then is scaled by a constant 
to normalize the model. We use the Altenhoff et al. ( 1994 1 



compilation because it represents the largest stellar sample 
available. We later (§5.3) compare a Tau to more recent 
interferometric observations by Cohen et al. (20051. 



Robinson, Carpenter, & Brown ( 1998 1 give estimates 



of rie in Q Tau for three epochs (1990, 1994, and 1996) 
based on near-uv C II spectra obtained with the Hubble 
Space Telescope: ric — 10^ cm"''. Unfortunately there were 
no contemporaneous lUE observations in 1986 when the 
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Table 1. Adopted properties of a Tau. 



Quantity 


Value 


Unit 


Spectral Type 




K5 III 




21.1 ± 0.2 


mas 


Teff 


3871 ± 48 


K 




44.56 ±2.94 


Rq 


logg* 


1.5 ±0.15 


cm s"'^ 




2.30 ± 0.85 


Mq 


[Fe/H] 


-0.15 ±0.20 






6.50 X IQ-i 


erg cm" 


pMg, Il/-Pbasal 


1.69 





Reference 

Mozurkewich et al. 
ditto 



(20031 



Decin et al. 
ditto 

inferred (g*) 
ditto 



(20031 



Martinez et al. 
ditto 



(2011 !+(/-. 



IRAM 250 GHz observations were made. Electron densities 
inferred from the C II] 2325 A multiplet ratios are weighted 
towards regions of high To and ne. Because semi-empirical 
models show that no is constant to within a factor of ~ 2, 
we adopt a mean that is half the peak observed value, 
< nl >= nl{HST)/A. 

The predicted flux is then obtained from equation (|3| 
with G = 1.08 and (TBr) given by the model. Combining the 
stellar parameters adopted in Table 1, we derive J250GHz = 
52mJy, likely a fortuitous agreement. The model predicts a 
radial optical depth of Ar,i=i = 0.05 and Ar^=i/2 = 0.09. 

4.2 Star to star variations in Ar 

Q Tau's Mg II h & k flux is about 1.7 times basal, and it 
might be slightly metal poor (see Tables 2 and 1, respec- 
tively) both of which affect the optical thickness. To scale 
this chromospheric thickness to other red giants, and remove 
explicit dependence on a Tau's stellar parameters, we define 
an optical thickness Atj^ (/i) for a fictitious star with a Tau's 
Tcff but now with basal-fiux heating [T =1) and solar met- 
alicity. Using the Ayres| ( |1979[ ) scaling of equation ([5| gives 

Combining equation (|5| and equation Q, the chromo- 
spheric optical thickness becomes 



reported in the two stellar samples, so we do not consider 
these objects further here. Because we are interested in the 
chromospheric emission, we need to avoid stars with sig- 
nificant ionized wind opacity, as is likely to occur at lower 
gravity (higher luminosity), e.g., the eclipsing binary C, Aur 
(K4 lb + B5 V) ( [Harper et "al]|2005[) and its spectral-type 
proxy A Vel (K4 lb) ( Carpenter et al. 1999 1, and earlier spec- 



0.04 Afo F-^ 



^^"^^^ M ApeQ i^basai V3871 Ky V 250 GHz , 
Note that the surface gravity term has now cancelled. 



V3871 Ky V 250 GHz y 



5 COMPARISON WITH OBSERVATIONS AND 
PREDICTIONS 

To test our semi-analytic model for the mm thermal fluxes, 
we require a sample of single inactive mundane red gi- 
ants with well determined stellar parameters. The 250 GHz 



(1.2 mm) single 30 m dish IRAM survey of Altenhoff et al 
( |1994[ ) presents the largest sample of measurements for 
these stars. We combined this survey with the accurate (j)* 
and Teff sample from the Mark III optical interferometer 
I Mozurkewich et al.||2003 l. The stars in common are listed 
in Table 2. We also have added a Lib (M3 HI) from the 
smaller sample of Dehaes et al. (2011 1. This combination of 
data enables us to make initial tests of the model. In Fig- 
ure 2, we compare the observed {Tbt) (G=1.08), with the 
predicted basal-limit for 250 GHz. 

The spectral-types of interest are single G, K and M gi- 
ants. Unfortunately there were no single G giant detections 



tral types. We also exclude the hybrid chromosphere stars 
7 Aql (K3 II) and t Aur (K3 II) because the hybrid a TrA 



(K3 II) shows evidence for a warm ionized wind (Harper 



et al. 20051. We include a Boo (K2 HI) even though it is 



thought to have a 10,000 K partially ionize d wind, becau se 

-10 ■ 



the mass-loss rate M = 2 x 10"'" Moyr"^ |Drake|[l985) is 
small, although perhaps not entirely negligible. For the later 
spectral-type semi-regular M giants it is not clear when the 
energy and momentum deposited by small amplitude radial 
pulsations begins dominate the sources present in earlier 
spectral-type red giants and begin to fundamentally alter 
the chromospheric structure, see for example [Eaton et al.| 
([1990|). We thus include R Lyr (M5 H I), while noting that 
[Luttermoser, Johnson, &: Eaton[ ( [1994[ ) found it was not pos- 
sible to construct a time-independent semi-empirical model 
for g Her (M6 HI). 

For the sample, we have estimates and errors for (j), T^a 
and F. In Figure 3 we directly compare the model (with 
G = 1.08) to [Altenhoff et al.| ([l994|). The model is broadly 



.(20) consistent with the observed fluxes but because of the large 
uncertainties (^ 15%) of the IRAM measurements it is not 
clear whether there is additional intrinsic scatter that is not 
accounted for in the model. Combined Array for Research 
in Millimeter-wave Astronomy (CARMA) and ALMA ob- 
servations with flux signal-to-noise ratios of ~ 20 will be a 
defining test of our simple model of time-independent but 
inhomogeneous red giant chromospheres. 



5.1 Spectral Indices 

Spectral indices, a (where cx 1^°'), are an empirical mea- 
sure to help interpret multi-frequency radio observation. For 
a constant temperature source with an angular diameter in- 
dependent of frequency a = 2, refiecting the frequency fac- 
tor in the Planck function. Many q's have been obtained at 
quite different frequencies where the fiux can arise from dif- 
ferent atmospheric layers (e.g., at cm- wavelengths these can 
reflect, say, the wind acceleration at 6 cm but the chromo- 
sphere at 2 cm), which limits their utility. With the greater 
sensitivity of the Karl G. Jansky Very Large Array (JVLA) 
and ALMA, a finer sweep of frequencies can be made and the 
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Table 2. Stellar properties, 250 GHz radio fluxes, and chromospheric optical thicknesses. The 
stellar properties are from Mozurkewich et al. (2003), the mm fluxes, except a Lib, are from 
Altenhoff et al.(1994) and the Mg 11 fluxes are from Martinez et al. (2011). 



Star 


HR 


HD 




^ cff 


rr(T a\ 


(h 




Flux 


C T i t LLJj I 


^ 1 basal 


At 








Tvnp 




(K\ 


1 'yn /I G 1 
I / / lllj j 


( Tti n Q\ 






I'Mp- it h k~l 


yfM — 1/ zj 


/3 And 


337 


6860 


MO 111 


3763 


46 


13.75 


0.14 


25.0 


4.0 


2.2 


0.15 


Q Get 


911 


18884 


Ml. 5 111a 


3578 


53 


13.24 


0.26 


15.0 


3.0 


2.0 


0.10 


p Per 


921 


19058 


M4 11 


3281 


40 


16.55 


0.17 


28.0 


3.0 


1.7 


0.05 


a Tau 


1457 


29139 


K5 111 


3871 


48 


21.10 


0.21 


51.0 


6.0 


1.7 


0.09 


t Aur 


1577 


31398 


K3 11 


4086 


50 


7.50 


0.07 


13.0 


3.0 


1.7 


0.18 


rj Gem 


2216 


42995 


M3 111 


3462 


43 


11.79 


0.12 


20.0 


5.0 


3.0 


0.13 


fi Gem 


2286 


44478 


TV T T T T 1 

M3 lllab 


3483 
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0.15 
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0.06 
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3705 


80493 


K7 111 
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47 
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0.07 
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1.9 
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3748 


81797 
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4069 
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0.17 
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4301 
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4637 


62 


6.74 


0.10 
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2.0 


1.0 


0.22 


a Boo 


5340 


124897 


K2 111 


4226 


53 


21.37 


0.25 


78.0 


8.0 


2.7 


0.11 


/3 UMi 


5563 


131873 


K4 111 


3849 


47 


10.30 


0.10 


16.0 


4.0 


2.2 


0.12 


7 Dra 


6705 


164058 


K5 111 


4013 


52 


9.86 


0.13 


10.1 


1.3 


1.2 


0.12 


52 Lyr 


7139 


175588 


M4 11 


3330 


44 


11.53 


0.16 


13.0 


4.0 


1.9 


0.07 


RLyr 


7157 


175865 


M5 111 


3174 


41 


18.02 


0.22 


14.0 


4.0 


1.5 


0.04 


7 Aql 


7525 


186791 


K3 11 


4099 


50 


7.27 


0.07 
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2.2 


0.24 


/3Peg 


8775 


217906 


M2.5 111 
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42 


17.98 


0.18 


23.0 


5.0 


1.6 


0.07 
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133216 


M3.5 111 


3634 


110 


11.00 


0.05 


12.1 


2.0 


2.0 


0.11 
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Figure 3. Observed minus model 250 GHz brightness temperatures. The model has been calibrated against a Tau and the uncertainties 
in the Mg 11 fluxes have not been included in the errors shown. While the model is consistent with the existing data, new higher accuracy 
observations are required to assess the model's viability as a mm flux estimator. 



resulting spectral energy distributions should provide more 
powerful diagnostic information. 

Using the model above we see that the 100 GHz flux 
is predicted to come from the upper portion of the chromo- 
sphere: by 250 GHz the chromosphere is becoming transpar- 
ent; while at 900 GHz the chromosphere can be become quite 
transparent, revealing the cold molecular layers above the 



photosphere, as seen in the CO fundamental bands ( Wiede- 



mann et al.|1994 l. For the sample in Table 2, we find spec- 
tral indices using ALMA Band 3 (100 GHz) and Band 9 
(660 GHz) lie in the range ~ 1.60 ±0.15 indicating that the 
mean temperatures are declining with increasing frequency. 
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Figure 4. Dehaes et al. (2011) 1.2 mm MAMBO II: Observations minus model brightness temperatures plotted with same y-axis as used 
in Figure 3. The SIMBA fluxes for a Boo and a Get lie outside this range. The model is consistent with these higher S/N fluxes, but 
note Dehaes et al. discount the MAMBO II fluxes for a Boo and a Get owing to unsuitable weather conditions. 



5.2 Variability 

The mm optical depths scale with chromospheric heating, 
so that observed changes in major cooling channel fluxes in 
the uv can reasonably be expected to induce corresponding 
changes in At, i.e., 5 J- oc / (5Ar). Few quantitative studies 
exist of the temporal chromospheric flux variability of red 
giants, although both a Boo and a Tau have been well ob- 



served with lUE. For a Tau, Cuntz et al. ( 1996 1 find a peak 



to peak Mg II h & k flux variations of ±4% on day time 
scales, whereas [Robinson, Carpenter, fc Brown| ( |1998[ ) find 
a slightly larger ~ 15% from three HST observations over 



seven years. McClintock et al. ( 1978 1 found little evidence 
for variability in a Boo's chromospheric emissions above 
Copernicus 's instrumental noise; 25% over tens of hours and 



20% on year time scales. Ayres et al. ( 1995 1 reported a 6% 



rms variation in the Mg II h line observed by lUEawex nearly 
a decade, although a smaller variability on day time scales 
was potentially masked by a long-term ~ 10% secular trend. 
[Brown et al.| ( [2008j find long term Ca II H & K variations at 
the 10% level with a suggestion of a > 14 year cycle, a Boo 
has a significant Mg II h & k emission above the basal-flux 
and it is possible that the excess magnetic heating is re- 
sponsible for changes in the chromosphere and mm-radio 
emission (see also discussion in |Sennhauser fc Berdyugina 
[2011[ ). [Judge et aL] ( [1993[ ) have examined the chromospheric 
flux variations of M giants, p Per (M4 II, ±10%), R Lyr 
(M5 III, ±4%), and g Her (M6 III, ±10%). These changes 
do not follow the variations in V magnitude. Perturbing the 
Ar's in the analytic model by ±15% imposes a fractional 
change in (TBr) (or 250 GHz flux) of < 5%. 

On the basis of these studies we do not expect signif- 
icant changes in the mm and sub-mm radio fluxes if the 
chromospheric material is uniformly distributed, such as may 



result from a high filling factor of the warm chromosphere, 
and/or spatial-averaged acoustic heating. 



5.3 Comparison with Recent Observations 

Red giants are feeble thermal radio emitters and obtain- 
ing accurate and precise fluxes is difficult. Some more re- 
cent observations have been made of red giant s by [Cohen[ 
et al.|(2005) (interferometric) and Dehaes et al. (20111 (sin- 
gle dish bolometer arrays). [Cohen et al.[ ( |2005| l have re- 
ported BIMA mm-fluxes for a Tau and a Boo obtained at 
217.8 GHz (1.38mm) and 108.4 GHz (2.77mm), and com- 
pared againt predictions based on photospheric Radiative 



Equilibrium (RE) values. Cohen et al.^ [2005) reported ev- 
idence for excess mm-emission from a Boo, which is con- 
sistent with the observed presence chromospheric heating. 
Further, the BIMA fluxes imply little change in temperature 
with frequency, whereas the present model reflects the de- 
cline of temperature as the higher frequencies probe closer 
to the photosphere. The 217.8 GHz flux for a Tau (when 
scaled) remains significantly lower than the IRAM 250 GHz 
flux. A comparison of the BIMA, RE and present model is 
given in Table 3. 

Dehaes et al.l (|20TT| present 1.2 mm SIMBA and 



MAMBO II observations for a sample of red giants. We com- 
pare these measurements to the model in Table 4. In general 
the agreement is quite good, except for a Get (M2 HI). The 
SIMBA fiux of 58.3 ± 4.1 mjy is much greater than their 
MAMBO II value of 23.6 ± 2.5 mJy, which the authors dis- 
count owing to unstable weather conditions. The later, nev- 
ertheless, is much closer to the IRAM and model fluxes. The 



SIMBA flux for a Boo is close to the Cohen et al. ( 2005 \ 



value (when adjusted for the different frequencies), whereas 
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Table 4. Observed and predicted 1.2 mm fluxes (mjy) for the standard stars of |Dehaes et al!] 
| |2011| l. Note that fluxes given within '[1' were discounted because of weather conditions. 



MAMBO 11 SIMBA IRAM Predicted 



Star Dehaes et al. Altenhoff ct al. 

a Boo {K2 111) [20.8 ±3.5] 105.6 ± 6.9 78 ± 8 > 57 

/3 And (MO 111) 23.5 ±2.7 < 40.1 25 ± 4 23 

Q Get (M2 111) [23.6 ±2.5] 58.3 ±4.1 15 ± 3 19 

/3 Peg (M2.5 111) 29.5 ±3.2 29.2 ± 2.5 23 ± 5 32 

t Aur (K3 11) 5.5 ±1.6 - 13 ± 3 ^8 

^ Umi (K4 111) 12.2 ±1.6 - 16 ± 4 13 

7 Dra (K5 111) 10.1 ± 1.3 - < 12 12 

cr Lib (M3 111) 12.1 ±2.0 - - 14 



Table 3. Observed and predicted mm fluxes (mJy) for a Tau and 
a Boo. Observed values are from [Cohen et al.| ( |2005| l and their 
predictions are based on radiative equilbrium models. 





a Boo (K2 111) 


a Tau (K5 111) 




217.8 


GHz 


Observed 


83.50 ± 17.1 


25.78 ±5.64 


Cohen predicted 


33.43 


32.38 


Present Model 


> 52.8 


40.8 




108.4 GHz 


Observed 


20.09 ±0.69 


13.97 ± 1.46 


Cohen Predicted 


7.94 


7.80 


Present Model 


> 17.7 


15.9 



our prediction is closer to IRAM. These measurements were 
made at different epochs, however. 

The limited data hint at variability at mm-wavelengths 
that is larger than one might anticipate based on the small 
observed chromospheric variability (peak-to-peak 15% that 
implies 4% in mm fluxes) 



6 DISCUSSION 

The time-independent inhomogeneous model appears to be 
consistent with most existing mm-radio observations. How- 
ever, because most measurements were made at a single 
frequency and have sizeable flux uncertainties, these data 
are only just beginning to provide insight into the chro- 
mospheres of inactive K and M giants. The characteristic 
temperatures of stellar chromospheres: Tj^^^, Tco, Tint, and 
Tchrom all are within 50% of the mean, and the unfortu- 
nate reality is that there is limited thermal dynamic range 
in this problem. The errors on (pt and Toff are at the few 
percent level so to test the limits of this model will require 
multi- frequency observations with S/N ~ 20 and as good 
an absolute calibration as possible (~ 5%), both of which 
are challenging. Multiple observations at closely spaced fre- 
quencies are highly desirable, because each flux would have 
a contribution from overlapping spatial regions and provide 
a much better measure of the mean temperature gradient 
than spectral indices based on widely separated frequencies. 
On this point we note that using Teft and (p, to estimate the 
mm-radio flux in equation (j4]| can give reasonable values 
even though it is a physically unsound model. 

We have not explored the G giants, but the model could 



be extended to include a filling factor of hotter material not 
at the boundary layer temperature ( Wiedemann et al.|1994 l 
by setting, for example, Tco = for that component. 



6.1 Variability 

Millimetre observations, such as those from JVLA Q-band, 
CARMA, and ALMA, have the potential to distinguish be- 
tween the debated acoustic or magnetic origin of chromo- 
spheric basal-flux heating. Numerical simulations show that 
granule sizes for main-sequence and giant stars scale with 
the surface pressure scale-height ( Freytag et al.|200"2 Nord- 
lund fc Dravins|[l990 l and the grain size is ~ lOHpres- The 
number of granules on the stellar surface will scale as 



: 4 X 10' 



6 / Rq \ 



cff.Q 



(21) 



which for a Tau gives ~ 10, 000 granules on the visible 
disk. This scaling appears consistent with findings from Ke- 
pler ( Mathur et al.pOll i. If the horizontal coherence scale 
of the acoustic shock heating is of the order of the gran- 
ule size ( [Loukitcheva et al.| 12008] ) then the spatially un- 
resolved radio fiux will have contributions from ~ 10, 000 
phase-uncorrelated patches. In the acoustic shocked heated 
interpretation, at any one moment there are likely to be sev- 
eral shocks passing through each patch, so one might expect 
that a basal-fiux red giant to show quite small variability 
at mm wavelengths in this scenario. Barring limb projection 
effects, the rms fractional fluctuations of the total intensity 
would be of order 1/-^10, 000 = 0.01. However, if observa- 
tions show significant variability, then it is likely that there 
is an additional, spatially more correlated source of chro- 
mospheric heating, i.e., magnetic active regions or possibly 
non-radial pulsations. To explore this aspect, it would be 
important to observe giant stars whose Mg II h & k or Ca II 
H & K fluxes are as close to the basal limit as possible. 



6.2 Molecular Component 

In the Sun the temperature minimum appears to be set by 
the balance between mechanical heating and cooling by H~ . 
Above this temperature minimum there appears to be an un- 
known fraction, perhaps as much as 80%, of the surface cov- 



ered by CO clouds or columns (Ayres 2002 I. In the red giants 
the CO, which appears to be close to radiative equilibrium 
temperatures, appears to cover most of the stellar surfaces. 



© 2012 RAS, MNRAS 000,[T|fT3l 



Chromospheric Thermal Continuum Millimetre Emission from non- dusty K and M Red Giants 11 



g Her (30 Her, M6 III, 



19 mas, Tefi = 3010 K) was not 



detected in the Altenhoff et al. ( 19941 survey (< 6mJy) and 
the upper-Umit to its brightness temperature is only 465 K. 
The present model predicts 28mJy and this is very dis- 



crepant. The star has been observed with both lUE (Judge 
et al.|1993[ ) and HST ( [Luttermoser et al.|1994| ) and it shows 
uv chromospheric emission, g Her and R Lyr both have sim- 
ilar and Tefi and have similar 2 cm radio VLA fluxes 
( Drake et al.|199T |. In the model presented here we adopted 
Too = CGTefi, but if the low fluxes for g Her and R Lyr are 
confirmed then the molecular component of the outer atmo- 
sphere may be even cooler, perhaps dominated by an H2O 
molsphere; and M5 III may mark the onset of a cooler domi- 
nant component. However, even this would not seem enough 
to explain g Her's low 1.2 mm flux and multi-frequency ob- 
servations are urgently needed. 7 Cru (M3 III) has been 
observed with the ATCA at 94 GHz . It has a large angular 
diameter {(j), ~ 2 5 mas, T^a = 3626 K; [Ireiand et al.|[2004 
Cohen et al.|1996 l and the present model predicts ~ 16mJy 



but the observed flux is significantly less (M. Cohen, priv 



It would seem that K and early M giants at the basal- 
Umit might be suitable mm flux calibrators at the 5% 
level. Random fluctuations in the brightness resulting from 
changes at the stellar photospheric granule scale should av- 
erage down, and chromospheric heating fluctuations derived 
from optical and uv emission suggest small variations in the 
mm-fluxes. For later spectral-types, the uncertain influence 
of this variable chromospheric heating is further reduced. 
That being said, the emission from the molecular material 
mixed in with warmer chromospheric plasma needs to be 
calibrated in some way, because it is challenging to estab- 
lish directly from theoretical models or simulations of stellar 
atmospheres. 

The semi-analytic model developed here is a time- 
independent 1-D model with all the inherent assumptions. 
When the next generation of theoretical chromospheric heat- 
ing models becomes available, one could analyse the thermo- 
dynamic structure and develop a new analytic representa- 
tion. Until then, the present model, based on the interpre- 
tation of empirical evidence should be regarded as a start- 
ing point. Even if real red giant chromospheres are highly 
time variable, the semi-empirical models upon which our 
mm radio predictor is based, still represent a measure of the 
temperatures, electron densities, and ionization required to 
produce the uv and optical spectral diagnostics. One could 
imagine an actual atmosphere as being a highly dynamic 
jumbled-up 1-D model, in which case the temperature rise in 
the warm chromosphere no longer is monotonic. The impli- 
cation for the radio fluxes would be to reduce them slightly 
as cool material is now present at greater altitudes. This 
systematic bias is being investigated. 

A test of chromospheric models for the largest angular 
diameter giants may be possible with ALMA, at frequen- 
cies above 250 GHz, in its most extended (16 km maximum 
baseline) conflguration. The chromospheres will be partially 
resolved but they will provide constraints on the millimetre 
limb brightening or darkening. 



7 CONCLUSIONS 

We present a semi-analytic 'toy' chromospheric model to 
examine millimetre thermal continuum radio fluxes of non- 
dusty and non-pulsating red giants. Equation ( |15| l shows, in 
a clear fashion, the contributions to the brightness temper- 
ature from the cool molecular layers at the base, and from 
the chromosphere itself, through weighting factors obtained 
from analytic solutions to the radiative transfer equation. At 
3 mm, the chromospheres are semi-opaque; while at 1.2 mm 
the chromosphere becomes optically thin and the cool molec- 
ular material becomes visible. The model predicts a hasal 
level of radio emission that depends on the temperature of 
the boundary layers. Studies at (sub-) mm will provide new 
insights into these especially enigmatic molecular regions. 
The model also permits an analysis of radio flux variability 
in response to changes in chromospheric heating as observed 
in Mg II h & k or Ca II H & K. If the chromospheric heating 
rate is not known for a given red giant, the model still can 
predict the basal mm radio fluxes. Significant variations in 
radio emission for basal- flux red giants, if detected, would 
indicate a source of structured large-scale heating, possibly 
magnetic in nature. The physical arrangement, e.g., columns 
and sheets, of the cool and warm chromospheric plasma can 
be constrained by the mm-radio limb-brightening of red gi- 
ants that can be spatially resolved when ALMA reaches 
its maximum baselines. A good test of this semi-analytic 
model would be to compare it against red giants with in- 
dependently determined T^s, 4>*, and Too and that can be 
monitored for changes in chromospheric heating. 

In summary, the present analytic model is consistent 
with most existing radio measurements and also highlights 
some possible anomalies. The sensitivity of mm-fluxes to 
changes in chromospheric heating can be evaluated and spec- 
tral indices are predicted for the existing thermal gradients. 
When a sample of multiwavelength (sub-) mm-radio obser- 
vations with good signal-to-noise becomes available these 
can be used to test the overall validity of the model. If the 
model passes this test, i.e., if there are not signiflcant star 
to star discrepancies, the observations can then be used to 
refine and test the mean thermal gradient in the stellar chro- 
mospheres. The spectral shapes could be used to better de- 
scribe the outward temperature gradient, or alternately they 
could invalidate the present model completely if there is no 
outward temperature increase. More importantly, these im- 
portant new radio observations can test new 3-D theoretical 
heating models when they become available. 

There appear to be a few stars with significant 1.2 mm 
flux differences between widely separated epochs and differ- 
ent observing facilities that appear at first sight to be larger 
than one would expect based on observed changes in chro- 
mospheric heating. These would suggest significant magnetic 
contributions for stars at twice the basal uv-fiux levels. We 
also have identified a few of late M giants which, while at 
the limit of applicability of the model, seem to have very low 
observed fluxes. However, we must await accurate new ob- 
servations of these stars before we can draw any signiflcant 
conclusions. 

While optical angular diameters and effective tempera- 
tures for many giants are known to a precision and accuracy 
of a few percent, we have shown that the physical dynamic 
range of Tbi for the mm-radio problem is limited. Future 
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progress will require high S/N multi- wavelength (sub-) mm 
radio fluxes with as small absolute flux calibration errors as 
possible. 
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